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Composite power system reliability involves assessing the adequacy of 
generation and transmission system to meet the demand at major system load 
points. Contingency selection was being the most tedious step in reliability 
evaluation of large electric systems. Contingency in power system might be a 
possible event in future which was not predicted with certainty in earlier 


research. Therefore, uncertainty may be inevitable in power system 


operation. Deterministic indices may not guarantee the randomness in 
reliability assessment. In order to account for volatility in contingencies, a 
new performance index proposed in the current research. Proposed method 
fares assimilates the uncertainty in computational procedure. Reliability test 
Probabilistic performancelndex systems like Roy Billinton Test System-6 bus system and IEEE-24 bus 
Contingency analysis reliability test systems were used to test the effectiveness of a proposed 
RBTS method. 
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1. INTRODUCTION 

Electricity is one of the greatest innovations of mankind. With unprecedented progression in 
population worldwide it is indispensable commodity of modern days. The economic growth of a nation is 
strongly coupled with its per capita energy consumption, which is only possible with a robust power system. 
Reliability is the measure of robustness of electric power system. Reliability, Security, and Stability were 
closely related terms often misinterpreted among researchers. As commendably outlined in [1], Power system 
reliability completely depends on its security, which further depends on its stability. Hence, security and 
stability were subsets of reliability. Security and stability were assessed by verifying the functioning of the 
system for a given set of operating conditions. However, reliability was being judged only by considering the 
time averaged performance over a reasonable period of time, generally over a year. So, security and stability 
were time varying attributes whereas reliability was time-averaged characteristic of electric power system. 

Reliability assessment of bulk electric power system comprises the following steps: 
a. Selection of system states 
b. Evaluation of adequacy of states 
c. Computation of reliability indices 

State space is gigantic for a practical interconnected power system. For instance a power system 
with ‘n’ components has 2n states. It may be highly cumbersome task to select and assess the enormous state 
space, which necessitate extreme computational endeavor. Several practises were proposed in earlier research 
for condensing the state space. Contingencies were ranked based on the severity of voltage degradation to 
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condense state space [2]. Bound estimates for evaluating the extremity of the contingencies may further 
lessen state space [3]. Misclassification of state space was considered as the masking effect in contingency 
screening algorithms. Research carried out by [4] deploys an adaptive procedure for compensating the 
masking effect in contingency classification heuristics. A contingency ranking table, indicating more or less 
events, for minimizing state space was employed by [5]. Chan KW et.al, [6] developed an online algorithm 
of system instability studies for considerable amount of reduction in state space with dynamic security 
contingency screening. 

Singular value analysis of load flow jacobian was used for ranking the contingencies with regard to 
voltage collapse to reduce the state space effectively [7]. Idea of integrating operators experience in 
contingency selection for decreasing the state space was proposed [8]. Artificial neural networks were also 
explored for contingency ranking. Fuzzy logic based MW ranking for contingency selection was explored 
[9]. He J et.al, [10] took an advantage of combining state enumeration and Monte Carlo simulations, by 
exercising state space partitioning method. Further, feed forward neural network was employed for online 
contingency screening [11]. Support vector machines along with scalar performance indices were exploited 
to rank contingencies according to the severity [12]. Ali Mansouri et.al, [13]investigated the impact of direct 
load control on nodal reliability of a deregulated power system. Saeed Afshar et.al, excersised preventive 
maintenance for electrical distribution system deliberating ageing [14]. 

Earlier research was centered on computation of deterministic performance indices and development 
of ANN based manipulating techniques for those indices. But, contingency by definition has a great degree of 
uncertainty imbibed in it. Previous methods employed deterministic indices which does not guarantee 
randomness in computation, which was the major problem never been addressed. Proposed Probabilistic 
index or Probabilistic Performance Index (PPI) integrates randomness in the process of its computation. In 
the current research investigations PPI method is proposed, which may be practical for contingency screening 
of power system. 


2. COMPOSITE POWER SYSTEM RELIABILITY EVALUATION 

There were two conventional methods available for calculating the reliability of composite power 
system: 
a. Monte Carlo Simulation Techniques 
b. Analytical Methods 

Monte Carlo methods of simulation computed the reliability indices by considering random 
characteristic of the system. These methods considered the problem as a sequence of experiments. In 
analytical methods, a mathematical model was being used to characterize the system and reliability indices 
were computed from the model. In the current research paper, state space enumeration approach is adopted. 
Analytical methods and Monte Carlo methods were having own advantages and disadvantages as discussed 
in [15], [16]. 

Computational time required for obtaining indices was relatively lesser in analytical techniques. 
Particular simplifying assumptions may be needed when the system became complex and large. Also, 
simulation methods may need long time for the computation. 


2.1. State space enumeration method 

Set of all possible outcomes of a random phenomenon can be designated as state space. In other 
words, combination of individual states of various components like generators, transmission lines, and load 
may be termed as state of a power system. The state space of a power system is represented as vector X. 


X={x;, > STOT Xp} (1) 


where, X1, X2,..., Xn are states of individual components like generators and transmission lines 

A system state in state space enumeration method may be designated by present state and transition 
that can occur between states. System state indicates the present working condition of working, failed, 
maintenance, etc. Change in state may cause the transition from one state to another. Such possible 
combination of all states makeup the system state space. Ideally, for precise assessment of power system 
reliability, entire contingency states of all orders needs to be considered. But it might be practically 
impossible task for assessing all contingencies. For example, RBTS-6 bus test system with 11 generators and 
9 transmission lines, total number of states needed are 27° (i.e. 10,48,576) for the assessment. 

To minimize the effort of a tedious task, numerous approximations and methods were developed. 
Truncation of State-Space (TSS) was proposed in the earlier research to decrease the computational burden 
[17]. In the TSS method, states with insignificant probabilities were excluded in state enumeration procedure. 
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Further, truncation is carried out on supposition considering that the probabilities of states denoting a large 
number of overlapping failures and lower number of failures. It was contemplated that probabilities of states 
with overlapping failures may be negligible when compared with probabilities with lower number of failures 
[1]. It was observed a great reduction in the number of states to be assessed and possibility for assessment in 
large variety of systems with the aid of suitable computer program [17]. Earlier, conventional methods of 
contingency ranking used the concept of Performance Index (PI) as described in Equation 2. 


Baga 2n 
PL = iL, [=| 2) 


where, 

Priow =the MW flow on line ‘P 

P/"**=maximum flow on line °P 

Once power flows were under the limit (rated line capacity), then larger values of ‘n’ leads to 

smaller value for PI. PI may be greater when one or more lines are overloaded in the power system. 
Contingencies were ranked based on the magnitude of the PI values. Contingencies below the threshold were 
discarded due to insignificant contribution for reliability. However, the earlier concept of performance index 
suffers certain shortcomings. The real power values used for computation of PI were absolutely deterministic. 
There was no uncertainty in earlier computations. Also, the computation of PI was dependent on exponent 
selection. Smaller values of exponent caused the masking effect. Greater the value of exponent lesser was the 
chance of misclassification of states and that was considered as masking effect [4]. In addition, larger value 
of exponent gives rise to nonlinearity. Accordingly, former research investigations disclosed that optimum 
integer selection for exponent will highly influence the PI calculation and further affects the contingency 
ranking [18]. 


3. PROPOSED METHODOLOGY 

In the current research investigations, a new index is proposed that overcomes aforementioned 
shortcomings. Though the computational procedure is laborious vis-a-vis conventional Performance Index 
(PI), the ranking of contingencies is quite promising. The proposed index is designated as Probabilistic 
Performance Index (PPI) and is defined in equation 3. 


PPI = Shr cow (Q EET 3) 


j=l Y Lmax 


where: W =Real non-negative weighting factor 

Y _=Order of the exponent 

N =No. of buses 

n =No. of contingencies 

Lij =MW curtailed at k" bus for j" contingency 
F; | =Frequency of failure for j contingency 

Value of the weighting factor (W) may be taken as 1. Order of exponent (Y) has to be optimally 

selected. Lower or higher values for Y may result in masking and nonlinearity problems in ranking the 
contingencies, respectively. The procedure for computation of PPI is illustrated in the following section. It 
was proposed to carry out research studies in a systematic manner as shown in Figure 1. Initially it was 
planned to gather the reliability data, component failure rate, repair rate for the test power system network. 
Henceforth, generator and transmission availabilities were calculated [1]. Further, state probabilities and state 
frequencies were computed and available capacity (Cav) was calculated for different states. Finally, PPI is 
determined and contingencies were screened for composite power system. 


4. COMPUTATION OF PROBABILISTIC PERFORMANCE INDEX 

In this section, the process of PPI calculation is illustrated by considering a small 3 bus Roy 
Billinton Test System (RBTS) [1]. Table 1 shows summary of generation data and Table 2 shows summary 
of transmission data for Roy Billinton Test System (RBTS). 
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Table 1. Generation Data for RBTS Example Table 2. Transmission Line Data for RBTS Example 
No. Capacity N ik Connected ji 
Plant eae (MW) (Ey) (hyo) Line an : (Ey) (yr) R(Q) XQ) B/2 
1 4 20 1 99 1 1 2 4 8 0.0912 0.48 0.0282 
2 2 30 3 57 2 1 3 5 8 0.08 0.5 0.0212 
ma sa 9° n 9 = 3 2 3 3 10 0.0798 0.42 0.0275 


Read System Reliability Data, 
Component Failure Rate, Repair Rate 


Calculate Generator, Transmission line availabilities using 
r4 
A _ A 
a+ 44 a+ ge 
Compute State Probabilities P1, P2,..., Pn 
Compute State Frequencies F1, F2,..., Fn 
v 


Calculate Cav for all states 


No 
Yes 
Vv 
Report Loss of Zero Load 
Load Lkj Curtailment 
v 


( Update Lkj Values [e— 


v 


Calculate Probabilistic Performance Index for all states 


Update PPI 
Values for all 
contigencies 


Figure 1. Procedure for carrying out of proposed research 


4.1. State Values for Generators 

From the proposed generator reliability data of Billinton and Allan [1] as encapsulated in Table 1, 
the availability and unavailability values are calculated equations 4 and 5. Based on the availability and 
unavailability values for generators, state values are determined as shown in Table 3. 


A =A, =4 =A, A 99_ 0.99 (4) 
Atu, 1+99 


U, =U, =U, =U, A = = 0.01 (5) 
Atu, 1+99 
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Table 3. State values for generators 
State values for plant 1 (4 generators) 


Units Units Probability 
Down Considered 
0 Aj Av A3A4 0.96059601 
1 4*(U1A2A3A4) 0.03881196 
2 4*(U1U2A3A4) 0.00058806 
3 4*(A1U2U;3U4) 0.00000396 
4 U,U2U3U4 0.00000001 
State values for plant 2 (2 generators) 
Units Units Probability 
Down Considered 
0 Aji Ad 0.9025 
1 2*(U1A2) 0.095 
2 U,U, 0.0025 


4.2. State Values for Transmission Lines 
From the proposed transmission line reliability data of Billinton and Allan [1] as encapsulated in 
Table 2, the availability and unavailability values are calculated equations 6 and 7. 


A, fA 1093 0.99636033 (6) 
A +u, 4=+1095 
U,=— 2d Í  —0.00363967 (7) 


1 


A +u, 441095 


Similarly, A, U, As, and U3 are computed. Summary of availability and unavailability values for 
transmission lines are tabulated in Table 4. 


Table 4. Summary of Transmission Line Availabilities 


Line No. Availability Unavailability 
1 A,=0.99636033 0.00363967 
2 Ao=0.99545455 0.00454545 
3 A;=0.99658703 0.00341297 


4.3. State Probabilities 

Most of the probabilistic applications in power system reliability assessment assume that component 
failures within a fixed atmosphere are independent events. When component failures are independent, then 
system state probabilities may be available by the product of unit state probabilities. From the state values of 
generators and transmission lines, state probabilities are computed as follows. 
State 1: with all components in service 


P,=(0.96059601x0.9025x0.98844633)=0.85692158 (8) 
State 2: with G1 out, all lines in service 
P,=(0.03881196x0.9025x0.98844633)=0.03462309 (9) 


Similarly, other state probabilities are calculated and tabulated using availability and unavailability values for 
transmission lines as shown in Table 5. 


4.4. State Frequencies 

During steady state, frequency of encountering a state is equal to frequency of exiting from the state. 
This concept is known as frequency balance approach [19]. Using frequency balance approach, state 
frequencies may be computed as, 
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F;=P;* (transition rate) (10) 


Also, from the state values of generators and transmission lines, the state frequencies are computed as 
follows. State 1: with all components in service 


Fy=P, * (Angi +2 Agot Arrt not Az3)=0.85692158(4+6+4+5+3)=18.85227476 (11) 


State 2: with G1 out, all lines in service 


F =P; * (3Ac1+ u g1 +2 Acat Art Arzt Axa)=0.03462309 (3+99+64445+3)=4.15477080 (12) 


Similarly, other state frequencies are calculated and tabulated as given in Table 5. 


Table 5. State Values for Transmission Lines 


Lines Out Lines considered Probability 
0 A1A2A3 0.98844633 

1 UAA; 0.00361076 

2 AUA; 0.00451345 

3 A1A2U3 0.00339509 

1,2 U,U2A3 0.00001649 
1,3 U,A2U3 0.00001237 
2,3 A,U2U3 0.00001546 
1,2,3 U,U2U3 0.00000006 


Amount of load curtailed for each contingency may be obtained using capacity available for each 
contingency and maximum load curtailed equation 13. 


Lxj = Cav = Lmax (13) 
where: Le; =Load curtailed at ‘k" bus for ‘j™ contingency 
Cav =Capacity available for each contingency 
Lmax =Maximum load at bus ‘k’ 


5. SIMULATION RESULTS AND DISCUSSION 

From state frequencies and load curtailments, the Probabilistic Performance Index (PPI) is 
calculated using equation 3 and summary of consolidated state values is portrayed in Table 6. Table 6 shows 
the PPI values for various contingencies that are arranged in descending order. Similar results are perceived 
using MATLAB based simulation analysis. Contingencies are classified into three categories as Non Critical, 
Moderately Critical, and Most Critical as shown in Table 7. 


Table 6. Consolidated State Values for RBTS 3 Bus System 


State Elements State State Frequency Load Curtailed PPI 

J Out Probability (P) (Œ) (Lip 

1 - 0.85692158 18.85227476 0 0.000000 
2 Gı 0.03462309 4.15477080 0 0.000000 
3 Gi, Gi 0.00052449 0.11436062 15 0.000061 
4 Gı, G2 0.00364454 0.63414996 25 0.005193 
5 G,, Li 0.00012648 0.15329376 0 0.000000 
6 Gi, L2 0.00015810 0.19145910 29 0.000636 
7 G;, Ls 0.00011857 0.11774001 20 0.000115 
8 G2 0.09020227 6.85537252 5 0.024275 
9 Gz, G2 0.00237374 0.30858620 35 0.002410 
10 Gz, Ly 0.00032951 0.38783327 5 0.000078 
11 Gz, L2 0.00041148 0.48438029 29 0.004077 
12 Gz, Ls 0.00030891 0.29115559 20 0.000710 
13 Lı 0.00313030 3.48402390 0 0.000000 
14 Lı, Ly 0.00001430 0.03150290 55 0.000062 
15 Lı, Ls 0.00001072 0.02128992 35 0.000011 
16 Ly 0.00391288 4.35112256 29 0.328550 
17 Ly, L3 0.00001340 0.02659900 110 0.000177 
18 Ls 0.00293466 2.62652070 20 0.057013 
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Table 7. Contingency Grouping Criteria 


Red indicates the most critical contingency 
Yellow indicates moderately critical contingency 
Green indicates Non-critical contingency 


Based on the contingency grouping criteria, Contingencies are itemized in decreasing order of PPI values as 
shown in Table 8. Contingency No. 16 (L2 out) is identified as the most critical contingency and 
Contingency No. 13 (L1 out) is recognized as non-critical contingency. 


Table 8 Contingency Ranking Using PPI 


State Elements State State Frequency Load Curtailed PPI 

j Out Probability (P)) Œ) 

16 L2 0.0039129 4.3511226 
18 L3 0.0029347 2.6265207 
8 G2 0.0902023 6.8553725 
4 G1, G2 0.0036445 0.63415 
11 G2, L2 0.0004115 0.4843803 
9 G2, G2 0.0023737 0.3085862 
12 G2, L3 0.0003089 0.2911556 
6 G1, L2 0.0001581 0.1914591 
17 L2, L3 0.0000134 0.026599 
7 G1, L3 0.0001186 0.11774 
10 G2, L1 0.0003295 0.3878333 
14 L1, L2 0.0000143 0.0315029 
3 G1, G1 0.0005245 0.1143606 
15 L1, L3 0.0000107 0.0212899 
1 - 0.8569216 18.852275 
2 G1 0.0346231 4.1547708 
5 G1, L1 0.0001265 0.1532938 
13 L1 0.0031303 3.4840239 


Red indicates the most critical contingency 
Yellow indicates moderately critical contingency 
Green indicates Non-critical contingency 


In order to validate the ranking procedure and to understand the contingencies graphically, RBTS 
test system is simulated using power world simulator [20]. Power world simulator is a cooperative power 
system simulation package intended to simulate high voltage interconnected power system. Figure 2 shows 
the simulation of RBTS system with all components in service. There are no line overloads and loss of load 
as shown in Figure 2. 

From the results of power world simulator, it is quite evident that outage of L1 does not cause any 
overload and loss of load on bus as shown in Figure 3. PPI ranking grouped such outage as ‘Non-Critical’ 
contingency. Simulation of Contingency No.16 illustrates an overload on line L3 due to outage of line L2 as 
shown in Figure 4. According to PPI ranking, this is the ‘Most Critical’ contingency. Hence, the PPI ranking 
is convincing in the par of power world simulation results. 

As per the amount of load curtailed, other contingencies like (L2, L3), (L1, L2), (L1, L3) may look 
like more critical, but their probability of occurrence is low. Therefore, from the results it was clearly 
understood that ranking of contingencies based on deterministic data is not appropriate. Such ranking may 
give deceiving results. PPI ranking was observed to be more promising as it takes out the uncertainty data of 
individual components into account. Further, PPI ranking limits the values between 0 and 1, in which any 
contingency may fall in the range. Also, ranking the contingencies based on other data like load curtailment 
may require wide range of values which again increase deliberately with size of the system. PPI ranking 
gives comprehensive scale for any dimension of the system. Furthermore, PPI ranking of contingencies for 
other test systems were simulated using MATLAB as shown in Figure 5 and obtained results are tabulated as 
shown in Table 9. 
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Figure 4. Illustration of Line 2 outage (Most Critical) 
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Figure 5. IEEE-24 bus Reliability Test system 


From the results, it can be clearly understood that the increase in bus number causes enlargement of 
state space. Although, PPI ranking gave proper results with greater bus size as shown in Table 9. 


Table 9. Results on Various Test Systems 


Total No. of contingencies Non a5 Most 
Testsystem Sidied Critica — CHtcal Critical 
RBTS 3 Bus system 18 
RBTS 6 Bus system 31 
RTS-24 Bus Reliability Test System 3000 


Red indicates the most critical contingency 
Yellow indicates moderately critical contingency 
Green indicates Non-critical contingency 


6. CONCLUSIONS 

This paper proposes a probabilistic performance index for contingency screening of composite 
power system. The modus operandi of the proposed index lies in translating probability data of the individual 
components into an index which incorporates the inherent randomness of the system. Conventional methods 
of contingency analysis use deterministic data which is quite irrational in reliability studies. As the index 
range from 0 to 1 the results are quite meaningful. The proposed approach has been tested on RBTS-3 bus 
system to validate the simulation results with rigorous theoretical calculations. The technique is also tested on 
large RTS-24 bus reliability test system to substantiate the results. The simulation results depict the 
appropriateness of the approach for classifying the contingencies into Non-critical, Moderately Critical & 
Most critical contingencies. Probabilistic Performance Index ordering facilitates the power system planner to 
identify the vulnerable areas and assists in developing reinforcement methods. The results also endorsed that 
the method is adoptable for large systems. 
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